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Abstract. We present a correlation between the templates was found by (Kogut et al., 1996a; Kogut et al 



ACME |/SP94 CMB anisotropy data at 25 to 45 GHz with 
the IRAS /DIRBE data and the Haslam 408 MHz data. We 
find a marginal correlation between the dust and the Q- 
band CMB data but none between the CMB data and the 
Haslam map. While the amplitude of the correlation with 



1996b). Significant correlations were found at each DMR 



frequency, but with a spectral behaviour in better agree- 
ment with free-free emission than with vibrational dust. 
These results have been confirmed by different experi- 



ments in various par ts of the s ky: Saskatoon ( de Olivcira 



the du pt is larger than that expected from naive models Costa et al., 1997[), OVRO ([Lcitch et al., 1997|) and 



of dust emission it does not dominate the sky emission. 

Key words: cosmic microwave background - Cosmology: 
observations 



1. Introduction 

The study of Cosmic Microwave Background (CMB) 
anisotropics has recently proven to be a powerful tool 



19 GHz ( ge Olivcira-Costa et al. ,1998;). At high G alac- 
tic latitudes, however, Python V ( |Coble et al., 1999| ) did 
not see any correlation. The regions of the sky covered by 
these experiments are shown in Fig. |l[ As noted above, 
free-free emission should correlate with Ha maps. Unfor- 
tunately, correlations found between Ha maps and the 
CMB indicate that the Ha-traced emission is too small to 
explain all of the observed correlation between CMB data 
and 100 /im data. 



(Drainc and Lazarian, 1998a; Drainc and Lazarian 



>■ , for observational cosmology ( de Bernardis et al,, 20Qflt 



1998b) have suggested an alternate explanation for the 



[Lange et al., 2000| ; |Hanany et al., 2000| ; [Balbi et al., 2000| ). 

One of the most important aspects of CMB analyses is to 
check whether the anisotropies observed are due to CMB 
temperature fluctuations or to foreground contamination 
such as diffuse Galactic emission. 

Diffuse Galactic emission is dominated at high fre- 
quency (above « 100 GHz) by thermal emission from 
dust. High quality tracers of this emission are given by the 
IRAS/DIRBE 100 /zm maps a nd have been made a vail- 
able in a user friendly way by ( Schlegel et al., 199S ). At 
lower frequencies, Galactic emission is dominated by syn- 
chrotron and free-free radiation. Our best tracer of the 
synchrotron emission of the Galaxy is given by the Haslam 
408 M Hz map ( | Haslam et al . , 1981) , Free-free emission is 



correlation between CMB data and dust maps, namely 
that this emission may be the result of rotational dust 
emission from elongated grains. This emission could be 
compatible with the observed spectrum of the dust corre- 
lated emission in the microwave frequencies. 

The sum of dust components, both rotational and 
vibrational, should show a local minimum at roughly 
70 GHz, increase to a peak around 10 GHz, and drop off 
at lower frequencies. Such behaviour has been observed by 



(de Olivcira-Costa et ah, 1999) in correlating the Tencr- 
ife 10 and 15 GHz data with the IRAS/DIRBE maps. 
This has been put forth as evidence supporting spinning 
dust as the origin of the dust-correlated emission. A re- 



cent analysis of yet more Tencrife data ( |Mukherjee et al 



traced I by Ha emission, but maps of this emission are not 
yet publicly available. 

Characterizing the correlations between these different 
Galactic components and data taken in the microwave is 
important in order to understand the spectral behaviour 
and the origin of Galactic emission which may contam- 
inate CMB measurements. The first significant, high-|6| 
cross-correlation between the COBE/DMR maps and dust 



2000D reports, however, that this result should be consid- 
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ered with care, as it originates mainly from a small, high 
Galactic latitude region. They finally say that the data 
does not allow any conclusion concerning the origin of the 
dust correlated emission. 

In this article we present a similar analysis using the 
South Pole 1994 (SP94) data ( pundersen et al., 1995] ). 
As correlations in the North Celestial Pole region have re- 
cently been questioned, we would like to obtain another 
measure of the correlation in a different region of the sky. 
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Ka — band coadded SP94 data 



Saskatoon 
o OVRO' 




Fig. 1. Regions of the sky (in Galactic coordinates) cov- 
ered by the different experiments that have been corre- 
lated with dust. Saskatoon is the grey region around the 
North Celestial Pole. OVRO observed 36 small patches 
of the sky at 2 degrees from the Pole. The main part of 
the Python V coverage is the grey region in the center of 
the lower part of the map. Tenerife observed a 10 degrees 
width strip around the North Celestial Pole (grey ring). 
COBE and 19 GHz both covered the whole sky. South 
Pole is at the right of Python V. 



Channel 


Frequency (GHz) 


Beam FWHM (deg.) 


K a l 


27.25 


1.67 


K a 2 


29.75 


1.53 


K a 3 


32.25 


1.41 


K a 4 


34.75 


1.31 


Ql 


39.15 


1.17 


Q2 


41.45 


1.11 


Q3 


43.75 


1.05 



Table 1. Frequencies and beam 
Maximum for the different channels. 



Full-Width-at-Half- 



More precisely, as correlations were found at low Galac- 
tic latitude (North Celestial Pole) and no correlation was 
found at high Galactic latitude by Python V, we wanted 
to explore the correlations in the high latitude regions that 
will be important for MAP and Planck. SP94 was a CMB 
experiment that used the Advanced Cosmic Microwave 
Explorer (ACME) which looked at such a high-|6| region. 

2. Data 

2.1. Microwave data 

The SP94 data were taken in two different bands with 7 
different channels: the K a -band with four center frequen- 
cies (27.25, 29.75, 32.25 and 34.75 GHz) and the Q-band 
with three different frequencies (39.15, 41.45 and 43.75 
GHz). A detailed description of the instrument can be 
found in (Meinhold et al., 1993), while the detail of the 



1994 observations and data reduction is described in (Gun 



dcrsen et al., 1995|) . Frequencies and beam width for these 
channels are given in Table [|. 
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Fig. 2. The top panels shows the coadded SP94 data 
for the Q-band and for the K a -band separately. The 
coaddition was performed using, for each pixel, the full 
covariance matrix (using the method explained at the 
end of section ||). The third panel shows the simulated 
IRAS/DIRBE data (dust) and the bottom panel shows 
the simulated Haslam data (synchrotron). For both simu- 
lated signals, only the 27.25 GHz band is shown, the others 
are very similar. 

2.2. Infrared data 

In this analysis, we used two templates for each channel : 
— For the thermal emission from dust, we used the 



IRAS/DIRBE map provided by QSchlegel et al., 1998| ). 
The initial resolution of this map is 6.1 arcmin FWHM. 
For the synchrotron emission, we used the Haslam 



maps (Haslam et al., 1981) whose resolution is 51 ar- 
cmin FWHM. 

Both maps were convolved to the SP94 beams (for each 
channel separately) using the beams given in Table |l| and 
then differenced, simulating the chopping strategy of the 
SP94 experiment. The simulated signals are shown, along 
with the SP94 data in Fig. ||. We see from this figure that 
the correlations, if non-zero, will be rather small. 
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3. Method 

For each of the seven frequency bands we have N = 43 
data points, all at a declination of —62° and ranging from 
23° and 67° in right ascension (see Fig. [l]). These data 
points were obtained using a sinusoidal 1.5° chop with 
smooth, constant declination, constant veloci ty scans. We 
directl y follow in th is section the notation of ( dc Olivcira- 
Costa |et ai, 199S() . We assume that the data is a lin- 



ear combination of CMB anisotropics x l CMB and M fore- 
ground components (including an unknown offset and gra- 
dient) that are given by simulated observations of the tem- 
plate dust and synchrotron maps using the SP94 beam and 
scanning strategy. In a vector-like notation, one gets : 



y = Xa + x C mb 



n 



(1) 



X is a matrix with N rows and M columns. It contains 
the template maps observed using the SP94 beams and 
scanning strategy, a it the M- vector we are interested in, it 
contains the correlation of the SP94 data with the different 
foregrounds contained in X . The experimental noise is n. 

The noise and the CMB anisotropics are each assumed 
to be uncorrelated Gaussian variables with zero mean. 
They each have non-trivial covariance matrices and to- 
gether a total covariance matrix of: 



C = (yy T )-(y)(y T ) 

= (xcmbXcmb) + (nn T ) = C C mb + C„ 



(2) 
(3) 



As discussed in (Gundersen et al., 1995), the noise in 
a given pixel in given band (K a or Q) is only correlated 
with the noise in the same pixel in the different channels 
of this band (4 channels for the K a band and 3 channels 
for the Q band). The noise from the K a and Q band are 
uncorrelated. 

The covariance matrix of the CMB is obtained through 
the window function of the South Pole experiment 



that can be found in ( Gundersen ct al., 1995| ) 
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where £ is the multipole index, i and j denote pixel indices 
and k and I denote the different channels. The window 
function was calculated accounting for the different beam- 
widths for each channel, and the underlying CMB sky was 
assumed to be the same for all channels. The beam width 
of each channel can be found in Table § C x is the CMB 
anisotropy angular power spectrum. We used the best fit 
fi = 1 model from ( Jaffc ct al, 2000] ) obtained with the 
combined COBE/DMR, BOOMERANG and MAXIMA-1 
data sets: Q tot = 1, A = 0.7, fl b h 2 = 0.03, n c h 2 = 0.17, 
n s — 0.975, tc — 0. As a check, we confirmed that the 
results do not change much using a flat band power spec- 
trum with Qrms-ps = 20, 25 and 30 fiK. The differences 
were small. 



We are interested in measuring the correlation coef- 
ficients with the various templates, we therefore want to 
measure a considering n and xqm b as noise (but account- 
ing for chance alignment between CMB and the templates 
trough the CMB covariance matrix). We therefore con- 
struct the following \ 2 '■ 



X 



(y 



Xa) T C~ 



1 (y - Xa) 



(5) 



Minimising this \ 2 leads to the best estimate of a : 

a=[X T C- 1 XY 1 X T C- 1 y (6) 

We make the fit simultaneously for all channels in order 
to account for the channel-to-channel correlation matrix 
elements but the templates (dust and 408 MHz emission) 
are fitted separately. The vector y contains the concate- 
nated data for all channels (N data points per channel). 
It is therefore a (n c han x N) vector. The coadded Q-band 
and coadded K a -band data are shown in Fig. ||. The co- 
variance matrix is an (n c ha<n x N) x (n c han x N) matrix. 
The X matrix contains the template simulated signal (one 
for each channel) as well as a constant term and a gradi- 
ent term for each channel, in order to account for the fact 
that such terms were removed from the SP94 data. The 
matrix X is therefore a (n c han x N) x (n c han x 3) matrix. 
In a given column (corresponding to a given channel), X 
contains mostly zeros, except for the rows corresponding 
to that channel where it contains the template simulated 
signal. 

The best estimate of a is given through the \ 2 mini- 
mization along with its covariance matrix: 



s = [x T c- 1 x] 



(<) 



We obtain a correlation coefficient for each channel but 
those are highly correlated according to the covariance 
matrix S. The correlation coefficients cannot therefore be 
averaged in a straightforward manner. We estimate the 
average of the correlation coefficients via: 



X 



(a - alfY,- 1 (a - al) 



(8) 



where 1 is a vector whose elements are all equal to 1 . The 
best estimate of a is then obtained by: 



Total(S- x a) 
Tota^E- 1 ) 



(9) 



where the Total function returns the sum of all elements 
of a matrix or vector. This is equivalent to a weighted av- 
erage, though in this case we also account for correlations 
between the different measures. The error-bar on a is then: 



a- a = V(Total(S-i))- 



(10) 



We do this for the K a - and Q-bands separately, and for the 
combination of the two. The results are shown in Table ^[ 
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Equations |]]| and |l0| implicitly suppose that the dust- 
correlated component has a flat spectrum (it has a spec- 
tral index n = 0). As was mentioned before, results from 
other experiments tend to favor spectral indices close to 
n ~ — 2 similar to free-free emission. In order to account 
for a non-zero spectral index, we model the emission as : 



v 

ao 1 — 



(11) 



where vq is a reference frequency and not an additional 
degree of freedom. Equation || then becomes : 



Ka bond 
Q bond 

all 


20.67 +/- 36.83 ^K„„ 
68.14 +/- 42.52 uK„„. 
22.19 +/- 36.76 /j.K c „. 


MJy/sr 
MJy/sr 
MJy/sr 


~' A 



X 



(a — aom) S 1 (a 



a m) 



where m is a vector with 7 elements m. 



(S)" 



(12) 



This 



35 

Frequency (GHz) 



X is not linear if one lets n vary as a free parameter. We 
therefore minimize it numerically for a grid of 300 values 
of n between -5 and 5. 



4. Correlations 

We compute correlations between the SP94 and both 
dust and synchrotron templates shown in Fig. [| and Ta- 
ble |2[ Those results are also plotted as a function of 
the frequency in Fig. |3|. The results we obtain show a 
marginal correlation of the Q-band data with the dust 
IRAS /DIRBE template whereas there is no correlation in 
the K a -band. No significant positive correlation is found 
in either band with the Haslam template. The uncertainty 
on the combined Q-band and K a -band correlation coeffi- 
cients is not as much smaller than that of the Q-band or 
K a -band alone as one might naively expect. In fact, this 
is due to the contribution of the CMB covariance matrix. 
For each channel, the CMB covariance matrix is compa- 
rable to the noise one and therefore combining the data 
reduces the noise contribution but not the CMB one that 



Ko bond 
Q bond 



-8.03 +/- 15.08 (jK,„.K J a 
-17.53+/- 15.18 ^K C „.K"' □ 
-12.41 +/- 14.59 ^K C „.K"' 
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Frequency (GHz) 



Fig. 3. Correlations coefficients between the SP94 data 
and the templates (top: with IRAS/DIRBE and bottom: 
with Haslam 408 MHz). The combined correlation coef- 
ficient assume a spectral index n = for the correlated 
component. 



( dc Oliveira-Costa et al., 199S ) and considered as a pos- 
sible evidence for the presence of spinning dust because 



becomes dominant. 



In Table |2j we also show the result for the analysis 
allowing the spectral index to vary as a free parameter. 
Combining the K a channels did not lead to any realistic 
value for the spectral index (the best fit, n = 5.0, was 
at the edge of the range we searched). Combining the Q- 
band data leads to a best fit for the spectral index of 
n = 0.8 ± 2.0. Combining all K a and Q data together 
leads to n = 2.1 ± 1.0. It is there clear that positive values 
are preferred. This is not a surprise, however, as the Q- 
band data obviously correlate with the IRAS data better 



of the fall at 10 GHz. The upper measurements ( Mukher- 
jee et al., 200d|) were published more recently and miti- 



gated the enthusiasm for spinning dust. As the values in 
this region are under question now, we plotted a square 
around these points. OVRO measured one point in this re- 
gion at 14.5 GHz ( Leitch et al., 1997 ) that seems to favour 
upper values of the correlation. Saskatoon data are taken 
from ( ie Oliveira-Costa et al., 1997 ) and show a la detec- 
tion in the K a -band (not significant) and a 1.9<r detection 
in the Q-band. The 19 GHz (whole sky survey) data were 
taken from (ie Oliveira-Costa et al., 1998). The result 



from PythonV at 41 GHz was taken from (Coble et al 



than the K a does. 

We summarize in Fig. |4| the measurements performed 
up to now of the correlation with the vibrational dust 
traced by IRAS/DIRBE lOO^rn : SP94 corresponds to the 
present article, the Tenerife measurements have been pub- 



lished in two different articles (de Oliveira-Costa et al 



199E; Mukherjee et al., 2000) reporting slightly differ- 



ent results obtained with almost the same data. We plot 
both measurements, the lower values were obtained by 



1999) and is fully compatible with no detection (the best 
fit is negative and does not appear in the plot). The COBE 
data points are from (Kogut et al., 1996b). The results 
quoted in the article were obtained by fitting DIRBE 
140/im to the DMR data. We therefore corrected them 
by the average ratio DIRBE 140^m/DIRBE 100/im to 
have them in the same units as the other d ata points. Th e 
FIRS point at 167 GHz was obtained by ( |Ganga, 199^) . 
We also overplotted in Fig. |] the predicted spectrum (in 
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Fig. 4. CMB/IRAS correlation coefficients as measured 
by various experiments. We overplot the spectra for vi- 
brationnal dust, spinning dust and the sum of the two. 
The vibrationnal dust spectrum assumes a spectral index 
of 2.0 and was normalized to match averages of DIRBE at 
Galactic latitudes higher than 20 degrees. The spinning 
contribution is normalized following (Draine and Lazar- 
ian, 1998a). To guide the eye, we also plot the best fit of 
the models to the dta in dotted lines. We also added a 
free- free emission spectrum with a spectral index of -2.1 
and an arbitrary normalization. 

terms of ratio to IRAS/DIRBE 100/im) for vibrational 
dust in green (we assumed a spectral index of 2.0) nor- 
malised with IRAS/DIRBE 100/im, and for spinning dust 
in blue. We also added the spectrum for free-free emission 
in light blue (arbitrary normalization and spectral index 



2.1). For t he spinning dust, we follow ( Draine and Lazar- 



ian, 19 38a) for the mixing between Warm Ionised Medium 
(WIM), Warm Neutral Medium (WNM) and Cold Neu- 
tral Medium (CNM) models using a respective fraction of 
0.14, 0.43 and 0.43. The res pective normalization of thes e 
models was also taken from (Draine and Lazarian, 1998a). 
The sum of both contributions is shown in Fig. |4| in red. 
The width of the curves for the spinning dust model is due 
to the galactic latitude dependance of the optical depth of 
the spinning dust components. We used all the latitudes 
more than 20° from the Galactic equator. The normal- 
ization is that given by (Draine and Lazarian, 1998a). As 
clearly yields too little emission to match the data, we 
have also done a fit to the points and arbitrarily raised 
the model curves by this amount. This is indicated by the 
dotted curve. 



5. Discussion and conclusion 

The signal to noise for the Q-band detection is only 1.6 
(89% confidence level) so this would not be considered a 
detection as most CMB experiments require at least 2a. 
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All 
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All 
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Table 2. Correlations coefficients between the SP94 data 
and the templates : We also give the coadded correlation 
coefficients found by combining the K a channels, the Q 
band channels and for all channels together. The second 
column gives the assumed correlated signal spectral index 
in each case (we considered n ^ only for the correla- 
tion with dust). The third column shows the coefficient 
of the correlation with the dust IRAS/DIRBE simulated 
signal (in /^K.(MJy/sr) _1 ) for each spectral index and the 
fifth column shows this correlation for the Haslam 408 
MHz simulated signal (in /jK.K~ 1 ). The fourth and sixth 
columns indicate the significance in terms of number of 
sigma of the correlation measured. No best fit was found 
in the range — 5 < n < 5 for the combined K a channels. 
For the combined Q channels, the best fit was found for 
n = 0.8 and is shown in bold. Combining all the channels 
leads to a best fit at n = 2.1 also shown in bold. 



We would note however that, as opposed to normal CMB 
anisotropy analyses, here we are comparing the data to a 
template and it is difficult to imagine that systematic ef- 
fects or analysis errors that would cause a random correla- 
tion with lOO^rn dust emission. The correlation could arise 
from random alignment between the CMB anisotropies 
and the dust template. This, however, should be taken 
into account in our analysis if our covariance matrices for 
the data and the CMB are correctly estimated. If these co- 
variance matrices were underestimated, the error bars we 
compute on the correlation coefficient (with Eq. ^ would 
be too small and therefore the significance of our result 
would be overestimated. In order to check the validity of 
our error bars, we correlated the SP94 data with tern- 
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plate dust maps obtained by rotating the initial template 
maps around the Galactic poles and by inverting North 
and South. The Galaxy was either rotated and/or inverted 
in 36 different ways (10 degrees each) to make 36 different 
simulations. We found a zero average correlation and the 
RMS of the correlations normalized by the error bars (cal- 
culated with Eq. appears to be 1.1. This shows that the 
covariance matrices of the CMB and the data and there- 
fore our error bars are correctly estimated, confirming the 
significance of the correlation coefficients we obtain. 

We do, however, find the following points interes ting: 
As seen previously ( |dc Olivcira-Costa et al., 1997 ) the 
correlation between the Q-band and the 100/im emission 
is stronger than the correlation between the K a -band and 
the 100/im emission. If the correlation we found is to be be- 
lieved, the ratio of the RMS of the 100/rni template times 
the fitted correlation coefficient to the implied sky RMS 
is 0.38 (in the Q-band). This result indicates that roughly 
14% of the power seen on the sky by ACME/SP94 Q-band 
could be due to Galactic emission. The Ce could go down 
by 14% and the amplitude by 38%. This however does not 
apply to the K a -band. This result is in qualitative agree- 
ment with (Gundersen ct al., 1995) and ( Ganga et al. 
1997|)7|both of which found different spectral indices for 
the K a - and Q-band data, though again, with low statis- 
tical significance. 

We have also done the above analysis using as a tem- 
plate not the raw 100/mr data but rather the extrapola- 
tions recommended by ( Schlegel et al., 1998| ) (model 8) to 
500GHz and 40 GHz. The \ z for the fit went down by a 
very small amount. The differences in the \ 2 were not sig- 
nificant and do not change the conclusions drawn above. 
This is not surprising as the SP94 data covers only a small 
region and the SFD model was designed to cope with the 
whole sky. 

In conclusion, we have shown that there is a very 
marginal correlation between the SP94 Q-band data and 
vibrational dust tracers, while no such correlation exists 
with the 408 MHz map. The amplitude of the correlation 



in the Q-band is larger than that seen by ( de Oliveira- 
Costa §t al., 1997|) . 
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